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The nonequilibrium free energy theorems show how distributions of work along nonequilibrium 
paths are related to free energy differences between the equilibrium states at the end points of 
these paths. In this paper we develop a natural way of barostatting a system and give the first 
deterministic derivation of the Crooks and Jarzynski relations for these isothermal isobaric systems. 
^ , We illustrate these relations by applying them to molecular dynamics simulations of a model polymer 

undergoing stretching. 
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INTRODUCTION 



Recently a set of revolutionary statistical mechanical theorems have been derived. These theorems involve dis- 
tributions of work and dissipation along thermodynamically nonequilibrium paths. The theorems are remarkable in 
that they constitute some of the very few exact statistical mechanical relations for systems that are possibly far from 
Q I equilibrium. Close to equilibrium one of these relations (the Evans Searles Fluctuation Theorem) can be used to 
Q ' derive Green-Kubo relations for linear transport coefficients. A second remarkable feature of these theorems is that 
they enable us to derive thermodynamic statements about small systems. The thermodynamic limit is not required. 
The third remarkable feature of some of the theorems is that one can derive exact relations about equilibrium free 
, energy differences by analysing nonequilibrium thermodynamic path integrals. 
I • The first relationship was presented in 1993 by Evans et al. [l| where an expression concerning the probability 
' distribution of values of the dissipation in the steady state was given, and tested. This work motivated a number 

■ of papers in which various fiuctuation theorems were derived, the first of which were the Evans-Searles Transient 
' Fluctuation Theorem Q, and the Gallavotti-Cohen Fluctuation Theorem Later the Jarzynski Equality (JE) Q 

~~' (also known as the work relation or nonequilibrium work relation) and the Crooks Fluctuation Theorem (CFT) [§| (also 
known as the Crooks Identity or Crooks Fluctuation Relation) were developed, giving expressions for distributions of 
time-integrated work along nonequilibrium paths that dynamically connect two different equilibrium thermodynamic 

■ states, and the free energy difference between the states. In the present paper we shall be concerned with the last 
^ ] two such theorems. 

I ■ There are proofs in the literature for the Jarzynski Equality in the isothermal isobaric ensemble for stochastic 
■ O ' dynamics 0| and for deterministic dynamics thermostatted homogeneously with a synthetic thermostat We 
5 , present the first derivation of the Crooks and Jarzynski relations for time reversible deterministic, systems which have 
a state point controlled by an external constant temperature, constant pressure reservoir. In order to carry out this 
derivation we surround a natural system which obeys Newtonian mechanics, the so-called "system of interest", with a 
reservoir region of fixed temperature and which maintains a constant pressure on the system of interest. This is a new 
model, and for many experimental systems, (i.e. those held at fixed pressure and temperature) it provides the most 
realistic arrangement yet available for simulating a process occurring under these conditions because it separates the 
Ci I reservoir from the system of interest. 

The existing statistical mechanical proofs 0, S, 0, of the JE [1, [ll| are for the change in Helmholtz free 
energy for a system held at constant volume and in contact with a thermostat. When applied to an experiment. 
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held at constant volume, the pressure of the initial equilibrium phase may be significantly different from the final 
equilibrium phase. Frequently experiments are performed at constant pressure rather than constant volume. Therefore 
it is often more appropriate to use the isothermal isobaric ensemble and the Gibbs free energy [H, [3| ■ Stochastic 
derivations, based on the Markovian assumption, exist for this case 0, 0, [ll|, • However, most actual experimental 
systems do not satisfy Markovian assumptions as the basic equations of motion are not stochastic (isl . IT^ . Stochastic 
derivations based on Markovian assumptions when applied to deterministic systems, are only valid for systems that 
obey the Langevin equation with white noise only \V^. Hummer and Szabo [3, present a formal derivation using 
Feynman-Kac theory and point out that their formalism remains valid if Hamiltonian or Schrodinger operators are 
used. At equilibrium such operators do not generate an ergodic canonical distribution because states with different 
energies never mix. When work is performed on the system they do not provide a mechanism for the dissipation of 
heat, and therefore although they are a useful theoretical construct for obtaining the JE or CFT, the nonequilibrium 
dynamics is not representative of a system, on which work is performed, which is in contact with a large thermal 
reservoir. The invocation of the initial canonical distribution implicitly assumes contact with a thermal reservoir. 
The use of thermostats provide a mechanism by which the .JE and CFT can be rigorously obtained for deterministic 
systems, and provides a realistic model of natural experimental systems where the system typically relaxes to the 
same temperature and pressure as it was at initially. Following the approach in @| we will derive the CFT [l^ and 
then the JE, which follows trivially. 

In the system we consider, the reservoir region obeys time reversible Newtonian equations of motion that are 
augmented with unnatural thermostatting and barostatting terms, while the system of interest obeys the natural time 
reversible Newtonian equations. We argue that as the number of degrees of freedom in the reservoir becomes much 
larger than that of the system of interest, and as the physical separation of these reservoirs from the system of interest 
becomes larg er the system of interest looses knowledge of the details of how the thermostatting and barostatting is 
achieved [20||. The system of interest simply "sees" that it is surrounded by a large time reversible reservoir region that 
is in thermodynamic equilibrium at a known temperature and pressure. The resulting Crooks and Jarzynski relations 
refer only to the work done on the system of interest and only the temperature and pressure of the reservoir. Hence 
we expect that the derived results of this gedanken experiment will apply to naturally occurring systems. While 
thermostats which only act on particles which form a solid wall enclosing the fluid are widely used, (for example see 
[20I. [2ll. [2^). the development of a natural barostat is new. 

Here we develop equations of motion for a fluid maintained at an externally controlled temperature and pressure. 
In contrast with previous work on wall thermostatted systems [20, 21, 22], the particles are not identifled as 'wall 
particles' or 'fluid particles', but the system and reservoir (which could be fluid and wall) are distinguished by their 
location. In fact, although the scheme can be applied to systems within thermostatting walls as before, it does not 
require any solid walls. As in previous work, the system of interest is purely Hamiltonian. Thus we arrive at the 
possibility of a fluid, large enough to contain a system of interest, which is, by the physical principle of locality, 
identical to a fluid regulated by a large external heat reservoir. However for this newly constructed fluid the effective 
decoupling from the thermal reservoir is axiomatic. This allows all of the previously developed work, on thermostatted 
dynamics, to be applied directly to this system. Thus given the principle of locality we obtain rigorous derivations 
for linear and nonlinear response theory, the fluctuation theorem, Le Chatelier's principle. Green Kubo relations, the 
CFT, the JE and the second law of thermodynamics [l^. [2^. 

II. THE EQUATIONS OF MOTION 

We wish to set up a dynamical model of a molecular system in the isobaric isothermal ensemble. The system 
of interest, where the equations of motion are Newtonian, will be surrounded by a reservoir region that uses Nose 
Hoover feedback [13, H,!!^ to control the pressure and the temperature. This feedback mechanism is unnatural, but 
deterministic and time reversible. The number of particles in the full system (reservoir plus system of interest), N, is 
constant. This can be achieved by making the total system periodic or by surrounding it by boundaries that prevent 
escape of particles. In the results presented below, we use periodic boundary conditions with a cubic unit cell of 
length L{t) centred as shown in Figure [H We fix the total momentum of the full system to zero {ptot = J2iLi Pi — ^) 
and the unnormalized centre of mass R of the full system to a fixed position which streams with the cell dilation, 
R(t) = X^iLi ^i'iiit) — R-(O) L(t)/L{0), by using appropriate equations of motion. Obviously if the centre of mass is 
set at the origin it will remain there, R{t) — Oy t, and for convenience we will restrict ourselves to this case. In order 
to specify the thermostatting and natural regions, we define qi as the magnitude of the displacement of particle i from 
the centre of its unit cell which in the case of the central cell is fixed at the origin qo = fi^ . In mathematical form 
Qi — KQi ~ <io)modL\- When the distance of a particle from the centre of its unit cell, qi, is less than r^, the dynamics 
are Newtonian (i.e. natural), see Figure [TJ 
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Figure 1: The diagram represents the three regions in the simulation. The inner darker circle represents the Hamiltonian region 
which is not affected by the barostat or the thermostat and whose diameter is held fixed, the white region is where the function 
g{q,V) increases from zero to unity and the outer lighter grey region is where g{q,V) is set to unity. The minimum width 
separating the white region from it's minimum image must be at least as wide as the maximum cutoff radius for the interaction 
between the particles. The outer boundary of the white region may fluctuate with the volume. 



The equations of motion are 



= pi/ni + avg{qi, V)c[i - -fq 
Pi = Fi ~ avg{q.i,V)pi - aTgiqi,V)pi - 
V = DayV. 



(1) 



where D is the Cartesian dimension. For a unit cell that is cubic, the time dependence of the cell length L ~ \/V 
is given by L = ayL. The volume and the temperature are controlled through the multipliers ay fear- These 
multipliers only operate if a particle is in the region qi > ri,. The function 5(qi, V) is a position sensitive switch that 
smoothly turns the thermostat and barostat off in the system of interest and on in the reservoir region: 



9{c^^,v{t)) 
5(q.,ni)) 



1 



Qi - n 



{L - r,)/2 - n 
1, q^>{L~ r,)/2 
0, qi < n- 



n<q^<{L- rc)/2 



(2) 



In this equation Vc is some distance equal to or greater than the longest distance over which the atoms interact. 
This means any work done on the system due to the volume V changing is accounted for by the equations of motion. 
Of course our choice for g{<\i,V{t)) is not unique. The total momentum is constrained as stated above by setting 
1p = jf ~ {o^v + otT)g{sij, V)pj), and is set to zero, and the system's centre of mass is constrained by setting 



7q — a.vY^^=i{9{lj^^) ~ SjLi ™j ) cind for convenience set to zero, R{t) = V This gives Newtonian 

equations including the streaming produced by dilation of the box when g[qi, V) =Q (i.e. in the central region), and 
is consistent with the usual isobaric-isothermal equations of motion when g{qi, V) — 1^ {qi, V) (l^. [27|. 

The time dependencies of the multipliers are such that the equilibrium equations of motion preserve the isothermal 
isobaric distribution function with fixed centre of mass at the origin and zero total momentum: 



exp(-/3(go(r) +PoF)J(ptot)(5(R) 
S^dV J dTeM-PiHoir) + PoV))6iptot)S{R) 



(3) 
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Here T =(q, p), /3 = l/(fc_BTo) where Tq is the thermostat temperature, Hq is the Hamiltonian associated with the 
standard Newtonian equations of motion and Pq is the externally appHed pressure. We note that this approach can 
be used to produce equations of motion that preserve the isothermal isobaric distribution while only considering 
configurational variables for the feedback mechanism jH] . The details of how the distribution function is preserved 
are given below, resulting in the following time dependencies for the multipliers 
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^Ev..(q.^)-q. + l-^+^U (4) 



where tt^tv are arbitrary time constants and ^ = DNkeT X^jLi 9{^i^ ' ■ Usually the forces all 

sum to zero — and thus ^ = 0, however if an external conservative potential is present (the case of 

a stationary optical trap for example j2§|) then ^ ^ 0. Note that as usual, care needs to be taken with the use 
g{qi,V)Fi • in a periodic system, and the correct expression for this case is given in the footnote {301- The 
constraints on the system, •jq and 7p may be made arbitrarily small by increasing the number of particles in the 
system or alternatively by choosing the arbitrary time constants tt and Ty to be large while keeping the number of 
particles in the system fixed. The phase space compression factor, i.e. the divergence of the equations of motion, 
A = Vq • q + Vp • p + dV/dV + dar/dar + day /day, is 

N N 

A = -aril - l/Af)DE 5(q., V) + ay{l - 1/N) ^ V,5(q., V) ■ + ayD, (5) 

i=l i=l 

and the extended Hamiltonian is defined as 

He = Ho + ^o^pT^kBTa + y^v'^v^bTq. (6) 

The extended instantaneous enthalpy is then defined Ie = He + PqV. We are now in a position to consider the 
time-dependence of fo{T,V) using the dynamics. We expect an equilibrium distribution function that is consistent 
with Boltzmann's postulate of equal a priori probability Eq. [H A dynamical system in equilibrium has the property 
that the distribution function for the ensemble is preserved df{T, V, t)/dt = 0, i.e. it is the time independent solution 
to the Liouville equation. Using the equations of motion it can be shown that 

f = W (7) 

and using, 

2m . -r . vr/(r, K - - A/(r, v, «). ,8) 

it is straightforward to show that this condition is obeyed when the extended distribution function is of the form, 
/(r',t) cx exp(— /3/£;(r')) where V represents the extended phase space vector T' = {T,V,aT,ay). Thus we obtain 
the standard isothermal isobaric distribution function with the addition of the two extra independent Nose Hoover 
variables, 

f{T') ^ /o(r, V) exp (4r^ + a^r^)) . (9) 

Given the standard isothermal isobaric distribution function, /o(r, V), Eq. [31 and assuming the system undergoes 
a uniform dilation it is easy to prove (i.e. the virial theorem) that at equilibrium. 



N~l N 



\ DV{t) 



i—1 



and that the kinetic temperature is given by the equipartition theorem, 

1 ^ 

ksT = — — y^(p^■p^)■ (11) 

D{N - l)m f^^^^^ ' ^ ' 

So we expect these equations to be obeyed if the above is implemented in a computer simulation. 
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Eq. [To] 
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Eq. [n] 


S.E. 


0.5 


0.50002 


0.0003 


1 


1.0006 


0.001 


3.5 


3.5001 


0.0015 


4 


3.9970 


0.003 



Table I: The average values obtained for the pressure and the temperature along with the standard error for a system with 36 
particles in 2 dimensions. 



III. TESTING FOR EQUILIBRIUM 

It is perhaps most effective to test for equilibrium using a small periodic system in two Cartesian dimensions. To 
do this we will make r;, = , so all particles lie in the reservoir region at all times. So, 

5(q„y(i)) = i (^l-cos(^2^-^*-^^ , < < (L - r,)/2 (12) 

g{q^,V{t)) = 1, q,>{L-r,)/2, 

where L = V^l^ and periodic boundary conditions are used [30| . 

We set the multipliers av = ax = 4, with the Cartesian dimension D = 2, the number of particles = 36 and 
used the WCA potential [31| with the Lennard Jones a and e parameters set to unity and the potential cutoff radius 
coinciding with rc in Eq. [I2l Two simulation sets were computed, firstly at the pressure Pq = 0.5 and temperature 
To = 1, then at the pressure Pq = 3.5 and temperature Tq = 4. The equations Eq. [10] and Eq. [11] were used to 
obtain the average pressure and temperature. At the two state points averages were calculated for 10^ time steps 
with dt = 0.001 using a Ath order Runge-Kutta integrator [sl]- This was repeated 50 times and the standard error 
was calculated. The results are shown in Table [B The difference between the input parameters and the averages is 
consistent with the standard error for the four variables. This provides strong evidence that the analysis we have 
presented so far is indeed correct. 



IV. THE CROOKS FLUCTUATION THEOREM AND THE JARZYNSKI EQUALITY 

The CFT 0, [l^ gives the probability of observing an amount of work AW = B, in transforming a system from 
an initial equilibrium state Ie,i to a final state Ie,2 (where the system will eventually relax to a new equilibrium), 
relative to observing the opposite amount of work AW = —B for the reverse process starting from the equihbrium 
state Ie,2- The ratio between the likelihood of these two observations is given in terms of the difference in free energy 
between the two equihbrium states. The JE [1, [lH gives the change in free energy between two equihbrium states in 
terms of the ensemble average of the exponential of the amount of work it takes to do the transformation. For both 
theories the two different equilibrium states must have the same temperature. These theories stand out as important 
due to the fact that they give differences in equilibrium thermodynamic potentials from nonequilibrium data. They 
also remain valid arbitrarily far from equilibrium. 

The Gibbs free energy is given by, 

G = -fcBrolnA(Po,To), (13) 
where the partition function A is the normalisation constant for /o(r, V) (see Eq. [3]) and thus. 



G — —ksTo In 



dV J dT5{ptot)SiR)eM-PiHo{T) + Po 



V)] 



(14) 



A natural definition of the rate of work (here the system's temperature Tq is held fixed) done on the system can be 
obtained by considering the first law of thermodynamics. 



dwt 



diE{r',t) _ diE{r'{t),t) 



dt 



dt 



(15) 



where the time dependence in the instantaneous enthalpy Ie is due to a parametric change, i.e. /_E(r', t) — lE(r', A(t)) 
in the Hamiltonian for the system. This definition remains useful when the parametric change alters the pressure 
Pq. If the time dependence is solely due to changing the pressure Pq then Eq. [TSl gives Wt = V dPf){t)/dt. Of course 
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if Pq is parametrically changed the blanketing technique is no longer able to circumvent thermostatting artefacts. 
At time t = the system is in an initial equilibrium with Ie{T',0) = Ie,i{T'), the system then undergoes some 
transformation during the time interval < t < t and then for times t > t the system relaxes to a new equilibrium, 
lE(r\t > t) — lE.2(r'), which is reached at the sufficiently later time tf ^ r. We need to calculate the probability 
of observing a certain amount of work done on the system in transforming it from an initial equilibrium with Ie,i 
compared to the opposite amount of work when the system undergoes the same transformation in reverse from an 
initial equilibrium with Ie,2- We consider an ensemble of trajectories, the number of which approaches infinity. 
Following the procedure exploited in the derivation of the Evans Searles fluctuation theorem j2;[23|) the probability of 
observing a trajectory with selected values of a property, such as work, can be obtained by summing over the product 
of the distribution function and the phase volume SVr for trajectories that satisfy the selection criterion. Hence, the 
probability ratio between the forward and reverse process is, 

Pf{Wr = B) _ Sr(0)|«i^=S,(iS 

5Vr{T'mh{T'm ^ ^ 

Pr{wr - -B) ~ Er(0)h.=B,.B^^r(r'.^(r))/2(,r'^(r))]' ^ ^ 

where /i is the distribution function for state 1 and /2 for state 2, the vector F''^ represents the time reversal 
transformation of T' (which has a unit Jacobian) i.e. p — > —p, ay —ay and ar —ax, and 5V{T'{0)) is 
the hyper-volume element in extended phase space surrounding the trajectory bundle. The summation is over all 
trajectory bundles which result in the amount of work done by the complete transformation process being in the range 
Wr = B± dB/2. It is also assumed here, and below, that all T lie on the ptot — and qcm — hypersurface. For the 
particular case here, using the equations of motion Eq. [T]this becomes. 



= B) _ Er 

5Vr{T'meM-PlEAT'm] exp(/3Gi) 
Pr{wr = -B) - Er'(o)|.„.=B,dB'5Vt(r''^(r))exp[-/3(/B,2(r^^(r))]exp(/3G2)' ^ ^ 

The quantities Gi and G2 are calculated from Eq. [T4]with Ie^X') given by either /^^(r') or /£;,2(r'). The flrst law 
of thermodynamics manifests itself in Eq. [15] and thus the rate at which heat is exchanged with the reservoir is given 
by ksTA. The probability distribution, in the streaming representation, starting from an initial distribution at time 
t = may be obtained from the streaming Liouville equation [H, [s^ and is 

/(r'(t),i) = /(r'(o),o)cxp(^-^*dsA(r'(s))^ . (i8) 

If we consider a packet of trajectories initially contained by the inflnitesimal element of volume SV{T'{0)) we may 
later observe them contained by the inflnitesimal element of volume SV{T'{t)) . By construction, the number of 
trajectories, or ensemble members, in these two inflnitesimal elements is conserved and thus we have 

sv{r'{t))f{T'{t)) = <5F(r'(o))/(r'(o)) (19) 

and in turn 

SV{r'{t)) = SV{T'{Q)) exp ^ ds A(r'(s))^ (20) 

Now by the flrst law Eq. [I5]we have Ie{t) = Ie{0) + Wr + ksT ds A{T'{s)): substituting this into the denominator 
of Eq. [17] along with Eq. [20] and noting /£;(r''^) = IsiT') gives, 

Pf{wT = B) _ J2r'(0)\w^=B,dB 

<5F(r'(0))exp[-/3/Ba(r'(0))] 



PriWr = -B) Y.T'm^.=B,dB '5^(r'(0)) exp[-/3(/B,i (r'(0)) + Wr 

where AG = G2 — Gi. Thus we obtain the CFT in the isothermal isobaric ensemble, 

Pf{Wr = B 



PriWr = -B) 



= exp[-/3(AG-B)] (22) 



which gives the probability of observing the amount of work Wr = B done in the transformation process from initial 
equilibrium state Ie,i relative to the probability of observing the amount of work Wr = —B for the reverse process 
starting from an initial equilibrium state Ie.2- It is now trivial to integrate Eq. [22] 

/C30 1*00 
dBpf{wr = B)exp[-(3B] = exp[-/?AG] / dBpr{wr = -B) (23) 
-00 J —00 
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Figure 2: The solid line is the change in Gibbs free energy as a function of distance between the two traps kr calculated from a 
series of equilibrium simulations using Eq. [27l The diamonds are the work done in stretching the polymer at a finite rate and 
the circles are the results obtained from the JE, Eq. 1241 

and arrive at the JE 

(exp[-/3w^]) = exp[-/3AG]. (24) 

Thus by taking the ensemble average of the exponential of the work done (defined through Eq. \TE\\ in transforming 
the system we obtain the change in Gibbs free energy AG. We again point out that this derivation, which starts from 
the basic equations of motion, has been carried out for a system of interest governed by the fundamental Newtonian 
equations of motion. This is only possible because the equations of motion we have introduced here allow us to in 
effect decouple the Newtonian system of interest from the larger thermal reservoir and take advantage of the principle 
of locality. This proof is vaHd far from equihbrium as opposed to existing stochastic proofs for the isothermal isobaric 
ensemble f3] which are only valid for Markovian systems. 

The JE and CFT provide practical approaches for obtaining AG by measuring the work done along an ensemble of 
nonequilibrium paths that transform the system from one state to another. However since the derivation of Eq. [24] 
is rehant on the specification of the equilibrium distribution function, numerical verification of this equation is also 
an indication that the distribution function indicated in Eq. [9] is actually generated by the equilibrium equations of 
motion. 



V. DEMONSTRATION: STRETCHING A CRUDE POLYMER MODEL 

To form a crude model of a single polymer suspended in a solvent of Cartesian dimension D = 2 we use the finite 
extendable nonlinear elastic (FENE) chain model ^S^l- The two end particles of the polymer are held by a pair 
of harmonic wells Uhix) = ^ktx^ , i.e. as in optical traps [2^. We may then apply the JE to the case where the 
polymer is stretched due to the two traps being moved apart. We note that previous work has addressed how to 
determine the free energy difference as a function of the polymer separation rather than the trap separation [18]. 
All particles, solvent and polymer, interact by the same repulsive WCA potential [311 and have the same mass. In 
addition neighbouring polymer particles are bound together by the FENE potential |.34 | 

Ufene = -\Krl In [l - (Ar/ro)^] , (25) 

where the particles must be separated by a distance Ar which is less than tq where the potential diverges and K is 
the strength of the potential. The spatial derivative gives us the force 

Ffene = -^Ufene = -K^vj (1 - (Ar/ro)2) . (26) 

We set the FENE parameters to K = lOe and Tq = 1.5cr where e and a are the parameters from the WCA potential. 
The total number of particles is = 300, the number of polymer particles is Ap = 8, the trap strength is kt = 2e/(j'^ , 




8 



the temperature is T = e/ks were is Boltzmann's constant and the pressure is Pq = 2tT^/e. The time is reported in 
units of a^Jmje where m is the particle mass and lengths are reported in units of a. The optical traps where initially 
set at a separation distance of Itr — 2.5 and stretched to a distance of Itr = 6.5. With Itr — 2.5 the average volume 
was found to be ly) ~ 521.9 ± 0.2 and with Itr — 6.5, ly) — 522.3 ± 0.2, where the figure in the error estimate is one 
standard error. The traps were centred abound the Newtonian region in the simulation. In the spirit of the method 
originally introduced by Hoover and Ree js^ the average equihbrium length of the polymer /p, in the direction of the 
vector connecting the two traps, was computed for 31 separation distances between and including the shortest and 
longest length. This data was then used to compute the equilibrium change in Gibbs free energy using the equation 

dG _ \ 

^ ^T^tryj-p — kr) 1 (27) 

with the trapezoidal quadrature. For the JE the work for Eq. [JHobtained from Eq. [T5]is computed from the equation 

dwt _ 1 dk , , , . 

— 2 dt ^'^P'^ "^P'' ' \ ") 

where Xp^r is the x component of the position of the particle at the right hand end of the polymer and Xp,i that at 
the left hand end. Both traps are located on the x axis. 

A total of 5 X 10^ nonequilibrium trajectories were computed starting from equilibrium with kr — 2.5 and stretched 
at a constant rate until Itr — 6.5 over a duration of r = 100. The results may be seen in Fig. [2l The average work 
done in the nonequilibrium stretching of the polymer is clearly greater than the change in free energy indicating that 
energy is being irreversibly dissipated into the thermal reservoir. It can be seen, despite the irreversible work, that 
there is excellent agreement between the change in free energy calculated using equilibrium methods Eq. [27] and the 
change in free energy calculated using the JE Eqs. [24lfc[28l 



VI. CONCLUSIONS 



We have introduced equations of motion that preserve the isothermal isobaric ensemble's distribution function yet 
still feature a region which is governed by the natural Newtonian equations of motion. This natural region can be 
chosen to be arbitrarily large. As we have argued previously (20l . |22 | if the thermostatted, barostatted blanket is 
far enough removed (from the system of interest) its details can not possibly effect the physical observations made 
in the system of interest. This amounts to the assumption of locality, which is one of the most important and well 
established assumptions in physics. When a far from equilibrium process occurs in the system of interest the heat fiux 
which is dissipated into the outer blanket diminishes with distance: for D = 3 this will be proportional to the inverse 
square of the distance. If the blanket is far enough removed from the system of interest the thermostat and barostat 
will only act on a region that is in local equihbrium and it is known that this class of thermostatted dynamics does 
not introduce any artifacts when in local equilibrium (l^ . 

In the arrangement described here, the system of interest is in a spherical region. Other physical arrangements 
could be designed, for example it could be a slit shaped region. This arrangement would be more useful for studying 
some nonequilibrium steady state dynamics (e.g. Poiseuille flow). 

Using these equations of motion we have derived the CFT and the JE in the isothermal isobaric ensemble that will 
be applicable regardless of how far from equilibrium the system of interest is driven. The only physical assumption 
necessary for this is the assumption of locality. Existing derivations for the isothermal isobaric ensemble using 
Markovian stochastic dynamics [6], can only be Hnked to the fundamental microscopic equations, at best, in the near 
equilibrium linear response regime. One of the remarkable features of these theories is that they remain valid far from 
equilibrium. We now have a proofs of the JE and the CFT that remain valid far from equilibrium in the isothermal 
isobaric ensemble. It is this ensemble, with the Gibbs free energy as the thermodynamic potential, which is most 
relevant to physical experiments and processes. 
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